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ABSTRACT 

We present a detailed investigation of the impact of astrophysical processes on the shape and am- 
plitude of the kinetic SZ (kSZ) power spectrum from the post-reionization epoch. This is achieved 
by constructing a new model of the kSZ power spectrum which we calibrate to the results of hy- 
drodynamic simulations. By construction, our method accounts for all relevant density and velocity 
modes and so is unaffected by the limited box size of our simulations. We find that radiative cooling 
and star-formation can reduce the amplitude of the kSZ power spectrum by up to 33%, or 1/zK 2 at 
I = 3000. This is driven by a decrease in the mean gas density in groups and clusters due to the 
conversion of gas into stars. Variations in the redshifts at which helium reionization occurs can effect 
the amplitude by a similar fraction, while current constraints on cosmological parameters (namely 
ag) translate to a further ±15% uncertainty on the kSZ power spectrum. We demonstrate how the 
models presented in this work can be constrained - reducing the astrophysical uncertainty on the kSZ 
signal - by measuring the redshift dependence of the signal via kSZ tomography. Finally, we discuss 
how the results of this work can help constrain the duration of reionization via measurements of the 
kinetic SZ signal sourced by inhomogeneous (or patchy) reionization. 

Subject headings: cosmology: cosmic microwave background — cosmology: large-scale structure of 
universe — galaxies: clusters: general — intergalactic medium 



1. INTRODUCTION 

In the last few years, significant progress has been 
made towards making precision measurements of the 
temperature anisotropy in the cosmic microwave back- 
ground (CMB) at sub-degree scales. Groun d-based ex- 
perim ents such as the South Pole Telescope dRuhl et all 
2004) and the Atacama Cosmology Telescope (jKosowskvl 
2003) have achieved the necessary combination of sen- 
sitivity and angular resolution to measure the primary 
CMB temperature anisotropy po wer spectrum to the 
seventh acou s tic pe ak and beyond (|Dunklev et al.l [20111 
iKeisler et aTJ 1201 ID . The PlanckB satellite will mea- 
sure the CMB temperature power spectrum at cosmic- 
variance limited precision to angular scales below one 
tenth of a degree. 

On angular scales smaller than ~4 arcminutes (£ > 
2700), the CMB power spectrum is dominated by 'sec- 
ondary' anisotropies, temperature fluctuations that are 
generated by the interaction of CMB photons with large 
scale structure between the observer and the surface 
of last scattering. The principal contribution to the 
secondary anisotropy signals comes from the Sunyaev- 
Zel'dovich (SZ) effect, which can be broken down into 
two components- the 'thermal' and 'kinetic' effects. The 
former describes the transfer of thermal energy from free 
electrons in the hot intra-cluster medium to CMB pho- 
tons via inverse-Compton scattering. CMB photons re- 
ceive a boost in energy, distorting the Planckian form 
of its spectrum. This gives the thermal SZ (tSZ) effect 
a unique frequency dependence which can be utilized by 
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experiments to extract this signal from the primary CMB 
signal and foregrounds. The kinetic SZ (kSZ) effect is 
caused by the dopplcr shifting of CMB photons via scat- 
tering off clouds of electrons with a non-zero bulk velocity 
(along the line-of-sight) relative to the CMB rest frame. 
Unlike the tSZ effect, the kSZ has the same frequency 
dependence as the primary CMB. 

The thermal SZ effect has now bee n detected for both 
large numbers of indiv i dual c l usters (IVanderlinde et all 
20101: [Marriage et all 120111 iWilliamson et all 120111 : 
Planck Collaboration et al.l 120111: iMarrone etaTl 120111) 
and as a seco ndary anisotropy signal in the CMB 
power spectrum ([Lueker et alj|2010t iDunklev et al.|[20Tl1 : 
iShirokoff et al.ll2011f ). While the kSZ effect has not yet 
been detected in either case, it is likely be done so first 
in the power spectrum. CMB temperature fluctuations 
sourced by the kSZ effect are proportional to the prod- 
uct of the electron density and line-of-sight velocity. The 
lack of an electron temperature weighting means that 
the contribution of low-temperature gas is more signifi- 
cant than for the tSZ effect (which is proportional to the 
product of electron density and temperature). On the 
other hand, the kSZ signal from individual groups and 
clusters is weaker. 

The kSZ pow er spectrum is also sen s itive to the details 
of reio nization ( Gruzinov "fcHiilpM IKnox et al I fl998L 
iMcQuinn et~aT1l2005l : iZahn et al.ll2005t llliev et al.ll2007f l. 
In models of inhomogeneous or patchy reionization - in 
which different regions of the Universe are reionized at 
different times - bubbles of free electrons around UV- 
cmitting sources are embedded in an otherwise neutral 
medium. If these bubbles have a large-scale bulk ve- 
locity then they will impart a temperature anisotropy 
onto the CMB. In this work we refer to the kSZ effect 
from reionization as the 'patchy' signal, while that from 
epochs after reionization is complete is referred to as the 
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' post-reioniza t ion' s ignal. 

iZahn et al.l (j2005l) and iMcQuinn et~all (|2005l) demon- 
strated that, to hrst order, the magnitude of the 
kSZ power from patchy reionization is dependent on 
the duration of reionization. Hence, when combined 
with measurements of the optical depth to reioniza- 
tion from the primary CMB power spectrum, the red- 
shift range spanned by the epoch of reionization can 
be constrained. For example, if reionization started 
at z = 14 and ended at z = 6, it would generate 
roughly 3^K 2 of patchy kSZ power (at I = 3000), 
while the range 8 < z < 12 would generate 1.5/iK 2 
(|McQuinn et al.ll2005fl . While there are several methods 
for probing the redshift at which reionization ended (e.g. 



Oh & Furlanctto 2005; Fan ct al. 2006; Lidz ct al. 2006; 



Becker et al. 2007; Bolton fc Haehnclt 2007; Ouchi et alJ 



2010HMortlock et al.l201ll and references therein), there 
are currently no other means by which the duration of 
reionization can be measured, making the kSZ power 
spectrum a unique and exciting probe. 

The patchy kSZ signal is expected to have a slightly 
different angular shape to the post-reioniza tion signal, 
peak i ng at a larger angular scale (£ « 2000. [Zahn et al.l 
120051 1201 ID . corresponding to the characteristic bubble 
size during reionization. In principle, the two compo- 
nents of the kSZ signal can be separated by a precise 
measurement of the power spectrum that encompasses a 
wide range of angular scales. In practice, however, the 
primary CMB signal swamps that of the kSZ at £ < 3000, 
while extra-galactic foregrounds dominate at i > 5000. 
Hence, differentiating the post-reionization kSZ power 
from that sourced by patchy reionization via their angu- 
lar scale dependence is an intractable task. To measure 
the patchy component of this signal it is therefore vitally 
important to have a good theoretical understanding of 
the post-reionization contribution. 

One of the principal aims of this work is to investigate 
the theoretical uncertainty in the post-reionization kSZ 
power spectrum. Specifically, we aim to investigate the 
impact of astrophysical processes such as radiative cool- 
ing of gas, the formation of stars and galaxies and feed- 
back from supernovae. We construct an analytic model 
for the kSZ power spectrum which we calibrate to the 
results of hydrodynamic simulations. These simulations 
are run in both the non-radiative regime, and including 
radiative cooling and star-formation. We are thus able to 
compare the impact of these processes on the kSZ power 
spectrum. Our model also enables us to investigate the 
cosmological scaling of the kSZ power spectrum. 

Several previous studies have measured the kSZ power 
spectrum directly from cosmological simulations by gen- 
erating synthetic sky maps, projecting thro ugh the simu- 
lation box stacked oyer multiple timesteps ([White et al.1 
[20021: lHallman et ail [20091 iBattaglia et al.1 120101 ) . The 
principal drawback to this approach is that very large 
simulation box-sizes are necessary to adequately capture 
the large-scale velocity flows that contribute significantly 
to th e kSZ signal at small angular scales (|Zhang et al.1 
2004). As we shall demonstrate, simulations with box 
sizes significantly less than 1 Gpc/h, systematically and 
substantially underestimate the kSZ power spectrum. 
Furthermore, high spatial resolution is required to re- 
solve small-scale baryonic processes such as cooling and 



star- formation. 

In this work we adopt a hybrid approach. We use high- 
resolution simulations to capture the effect of non-linear 
structure formation on the gas density power spectrum, 
and use the results to improve our analytic calculation 
for the kSZ effect. Our method accounts for all relevant 
velocity and density modes, and circumvents the pro- 
hibitive requirement of high spatial resolution and ex- 
tremely large cosmological boxes. 

This paper is organized as follows: in Section[2]we pro- 
vide an overview of the kSZ effect, describe our model 
for the kSZ power spectrum, and discuss how we im- 
plement the modifications required to account for the 
astrophysical processes in our simulations. In Section [3] 
we describe the hydrodynamic simulations used to cal- 
ibrate our model. In Section [4] we present results from 
the simulations, focusing specifically on the gas density 
and momentum power spectra and comparing these with 
the predictions of our models. In Section[5]we discuss the 
kSZ power spectra predicted by our models, their scaling 
with cosmological parameters, and the means by which 
they could be distinguished observationally. We also in- 
vestigate the impact of helium reionization on the power 
spectrum. Finally, we compare our model with the re- 
sults of previous work and the latest observations. 

Except when referring to specific simulations with spe- 
cific cosmological parameters, throughout this paper we 
assume a fiducial, spatially-flat, ACDM cosmological 
model consistent with the WMAP7 best-fit cosmological 
parameters, namely Hq = 71 km s — 1 Mpc _1 , £Im = 0.27, 
n b = 0.047, = 0.73, n s = 0.95 and cr 8 = 0.82. 

2. MODELING THE KINETIC SZ POWER SPECTRUM 

2.1. kSZ Basics 

Thomson scattering of CMB photons off clouds of free 
electrons with a coherent bulk velocity along the line- 
of-sight from the observer produce fluctuations in the 
observed brightness temperature of the CMB, 
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where <tt is the Thomson cross-section for an electron, 
x is the comoving distance to redshift z, n e is the free 
electron number density and v • n is the component of 
the peculiar velocity of the electrons along the line-of- 
sight. We are principally concerned with the kinetic SZ 
power spectrum in the post-reionization era, so the upper 
limit to the integral, z ro j, corresponds to the redshift at 
which reionization ends (i.e., hydrogen has been entirely 
ionized). Unless stated otherwise we assume a fiducial 
value of z ie i = 10. 
t is the Thomson optical depth, 



t(z) = <7 T 



n e {z') dx 
l + z'dz 1 



dz' 



(2) 



(3) 



where n e is the mean free-electron density, 

- _ XPg(z) 
fi e m p 

Pg(z) = /Og,o(l + z) 3 is the mean gas density density of the 
Universe at redshift z, and /i e m p is the mean mass per 
electron, where fi e = 1.14. We define x as tne fraction 
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of the total number of electrons that are ionized. We 
assume that, for z < z re ;, hydrogen is completely ionized, 
so x is dependent on the abundance and ionization state 
of helium, 

i-r„(i-JW4) (4) 

X l-F p /2 ' {) 

where Y p is the primordial helium abundance and N^ e 
the number of helium electrons ionized (which can be a 
function of redshift). Thus for Y p = 0.24, X ={0.86, 0.93, 
1} for neutral, singly, and fully ionized helium. For our 
fiducial model we assume \ = 0.86 (i.e., iV Hc = 0) at all 
redshifts. We explore the effect of helium reionization on 
the kSZ power spectrum in Section [SJ 

Writing n e = n e {l + S) we define the density weighted 
peculiar velocity q = v(l + 6), so that 
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KSZ temperature fluctuations are generated by the 
projected contribution of ionized gas with a non-zero pe- 
culiar velocity along the line of sight. A key property of 
the kSZ signal is that Fourier modes of q (q(k)) which 
have k parallel to n suffer severe cancellation when pro- 
jected along the line of sight. Therefore, only modes of 
q(k) parallel to n but perpendicu lar to k can contribute 
(see Uaffe fc Kamio nkowski 1998, for a rigorous demon- 
stration of this. Henceforth we use ~ to denote a Fourier 
space quantity). 

q can be decomposed into divergence free (qs) and 
curl-free (qg) components, which satisfy y.qs = an d 
V x q^ = 0, respectively. In the Fourier domain, q = 
<Le + <Lb, where q_E = k(q • k) and 



qs(k) = q-k(qk) 



(6) 



When projected along the line of sight, the peaks and 
troughs of the component of q parallel to k, qg, will 
cancel. Therefore, only the component of q p erpendic- 
ular to k, c{b, contributes to the kSZ signal (jVishniad 
[19871: Uaffe fc Kamionkowskl[T998l) . 

qs can be written as a convolution between the Fourier 
transform of the velocity and density fields, 
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where /i = k • k'. In the linear regime, v(k) is parallel to 
k, so only the v5 component of q can contribute to q^. 
The kSZ power spectrum is thus generated by the c ross- 
term 5v (|Ostriker fc Vishniadll986t lVishniadll987h . 

In the small angle limit, the kSZ angular power spec- 
trum can be calculated using Limber's approximation, 
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where k = t/x, A%(k, z) = k 3 P B (k,z)/(2-K 2 ), and 
the power spectrum of the curl component of the mo- 
mentum field, Pb, is defined by (27r) 3 PBi5(k — k') = 
<q B (k)q^(k')}. 



IVishniacI (119871) first calculate d an expression for A% 
(see alsolJaffe fc Kamionkow ski 1998; Dodclso n fc Jubasl 
fl995tlMafc Frvll2002l) . 
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where P55 and P vv are the linear theory density and 
velocity power spectra, and Pg v is the density-velocity 
cross spectrum (for clarity we have suppressed the red- 
shift dependence of A B (k) and P(k)). In the linear 
regime, the continuity equation allows us to relate the pe- 
culiar velocity field with density perturbations, v(k) = 
ik(/d/fc)(5(k), where / = dlogD/cHoga and D is the 
linear growth factor. Therefore, 



Plugging these into Equation [9j we obtain 
k 3 t d 3 V 
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where 
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Combining Equations [5] and H] gives the well-known 
Ostriker-Vishniac effect, which is the linear-theory part 
of the kinetic SZ power spectrum. 

2.2. Non-linear Contributions 

Several previous studies have investigated the impact 
of non- linear structure formation on the kSZ power spec- 
trum, showing that non-linea r corrections become large 
for angular s cales i > 10 00 Ex] [2 00f4 iMa fc Fry| [2002) 
IZhang et al.l l200l . IHtjI (I2000D and TMa fc Frvl (2002 s ) 



demonstrate that one can include the effect of non-linear 
density fluctuations to the kSZ power spectrum by ex- 
changing the linear theory matter power spectrum, Pgg, 
in Equation [TT] with the non-linear matter power spec- 
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IMa fc Frvl (|2002h argue that the kinetic SZ signal is 
less sensitive to non-linear velocity fluctuations than 
non-linear density fluctuations due to the 1/k 2 weight- 
ing in the former. Using hydrodynam ical simulations, 
iZhang et all (|2004h : iShao et al.l (|2011| ) show that, for 
k > 2 ftMpc^ 1 , the power in the curl component of 
the velocity field (generated by non-linear gravitational 
collapse) can exceed the linear theory prediction. To ac- 
count for this, they suggest an phenomenological correc- 
tion in which Pss(k') is replaced by its non-linear coun- 
terpart. In this work, we follow the non-linear correc- 
tion proposed by IMa fc' Frv ( 2002|) given in Equation [T3l 



4 



Shaw et al. 



However, we investigate to what degree velocity fluctua- 
tions in the non-linear regime may effect the kSZ power 
spectrum in Section [3] 

Throughout this work, we calculate the non-linear den- 
sity p ower spectrum, PH L using the HaloFit prescrip- 
tion of lSmith et"all ([20031) . We find that the dark matter 
power spectrum predicted by HaloFit is within 8% of 
that measu red in our non-ra diative simulations for k < 
4 /iMpc" 1 (|Rudd et al.ll2008D . although it systematically 
underestimates t he po wer spectrum at ft > 1 hMpc" . 
iHeitmann et al.l (|2010l) found similar results comparing 
HaloFit to a large suite of N-body simulations over a 
wide range of cosmological models. 

There is a much larger discrepancy between HaloFit 
and our simulation that includes cooling and star- 
formation (CSF). The formation of dense clumps of stars 
and gas at the center of halos produces a steepening 
of the dark matter density profile at small radii, in- 
creasing the density of the halo core. This process is 
sometimes referred to as 'ha lo (or adiabatic) c ontrac- 
tion' ([Blumenthal et al.lll986l iGnedin et alJl200l I2TJTT1 ). 
The matter power spectrum measured in dark matter- 
only simulations, or by halo models calibrated on such 
simulations (such as HaloFit), no longer matches that 
measured in simulations that includ e radiative cool- 
ing at wavenumb e rs ft > 1 ftMpc" 1 (iJing et al.l [20061: 
Rudd et al.1 120081 iDuffv et al.1 [20101 Ivan Daalen eTaLl 



We incorporate the effects of halo contraction us- 
ing the simple mo dification to HaloFit suggested by 
iRudd et"aH ([2008). This is implemented by multiply- 
ing the matter power spectrum by the ratio of Fourier- 
transformed NFW density profiles (Nav arro et al.l[T997t ) 
with two different concentrations, 

PuMcsl (ft, z) = $(fc, z) 2 PDMonly(ft, z) , (14) 



where 



$(fc,z) 



X(R v hk/c 2 ,c 2 ) 



A(i?virft/Cl, Cl) 



(15) 



A(i? v ir, k,c) is the Fourier transform of an NFW halo of 
virial radius i? v iJ3, and (02, c\) are the concentrations of 
the halo including and omitting the effects of radiative 
cooling and star-formation, respectively. We apply the 
correction to the matter power spectrum given by Equa- 
tion [TH using ci = 5, c 2 = 8.5 and R v i r (z = 0.55) = 
1.1 comoving Mpc/ft, (corresponding to a virial mass of 
M vir = 1.07 x 10 14 h- 1 M Q ). 

This procedure is clearly a very simplistic approxima- 
tion to the effects of halo contraction on the matter power 
spectrum. The correction to the power spectrum is deter- 
mined entirely by the effect of a change in concentration 
on the density profile of a halo with a characteristic mass 
and redshift. However, in Section 2] we demonstrate that 
it does provide a good agreement between our model and 
the CSF simulations at small scales. 

2.3. Thermal Pressure of Baryons 

4 The virial mass and radius are defined by Af v ; r = 
g7r.R*Jj r A c p c (z), where A c is the virial overdensity given in 

IBrvan fc Norman (1993) an d pc(z) is the critical density at red- 
shift z 

5 See IRudd et~aT1 l |2008h for an analytic expression for A. 
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Fig. 1. — The contribution, d? C 'i / 'dzdln(k) , to the total kinetic 
SZ power at £ = 3000 by density fluctuations at a given redshift and 
spatial scale. In each panel, the contours enclose the regions that 
contribute 25, 50 and 75 % of the total power (from inner to outer 
contours). The top panel shows mass and redshift dependence 
for the full kinetic SZ power spectrum including the non-linear 
density fluctuations. The lower panel shows the same contours for 
the Ostriker-Vishniac effect, in which only the linear density power 
spectrum is used. The dashed blue lines represent the size-distance 
relation, k = £/x(z), for I = 3000. Kinetic SZ power is generated by 
density fluctuations in the range 0.3 < k < 10 hMpc -1 , with non- 
linear fluctuations boosting the signal in the high-fc, low redshift 
regime. 

In previous work, it has commonly been assumed that 
<5gas = &m, that is, perturbations in the gas density 
exactly follow those of the dark matter at all scales 
(iDodelson & Jubasl Il995t Uaffe fc Kamionkowskl 119981: 
IMa fc Frvl I2002t IZhang et al. 2004). However, thermal 
pressure between baryons should erase density fluctua- 
tions in the gas distribution at small scales. Further- 
more, radiative cooling and star-formation can signifi- 
cantly modify the ga s density power spectrum at both 
large and small scales dJing et al.ll2006[ IRudd et alj|2008t 
iCasarini et al.ll201lUvan Daalen et al.ll201lD . 

In order to accurately predict the kSZ power spectrum 
we must be able to relate the power spectrum of gas 
density fluctuations, P gas (or more precisely, fluctuations 
in the number density of free electrons), with that of 
the dark matter. We incorporate the effects of baryon 
physics into a window function W(k), such that 
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W 2 (k,z)P^(k,z) 



(16) 



This then replaces the non-linear matter power spectrum 
in Equation [T31 Note that we assume the velocity field 
of the gas follows that of the dark matter, which is a 
reasonable assumption at large scales. 

In the non-radiative regime, the qualitative shape of 
W 2 (k) is simple to imagine. At large scales, and before 
the onset of gravitational collapse, we expect W 2 (k) w 1. 
However, at small scales, the gas thermal pressure force 
suppresses gas density perturbations and so W 2 {k) will 
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tend to zero as fc increases. Therefore, W 2 {k) acts as a 
filter, smoothing the gas density at some characteristic 
scale. As we shall demonstrate, the form of W 2 (k) for our 
simulation including radiative cooling and star-formation 
h as a more comp l ex dep endence on k. 

IGnedin fc Huil (|1998l ) demonstrate that, for coupled 
density perturbations in the linear regime, the ratio of 
gas to dark matter density fluctuations is well described 
by the form 

W G98 (k, *) = 0.5 (cxp (-k 2 /k 2 ) + [1 + 4(fe ) fc/)2]1/4 ) , 

(17) 

where the (redshift-dependent) characteristic filter scale 
kf is given by 

J__J_f t 2 , , D(t') + 2H(f)D(t') f* _dt^_ 
k}(t)~ D(t)J a{t)dt k 2 (t') J t ,a 2 (t") 

(18) 

Here fcj is the Jean's scale, 

kj(t) = -^-V^Gpm(t) , (19) 
cs{t) 

p m is the mean matter density of the Universe and eg = 
y/ dP/ dp is the mean sound speed at time t. As noted 
by IGnedin et al.l (|2003f ) . the filter scale at a given time is 
not directly proportional to the Jeans scale at the same 
time, but to the integral over the thermal history of the 
g as u p to that point. 

iHul (|200l and iMcQuinn etldl pQOl used Eq. Q7| to 
approximate the effects of thermal pressure in their cal- 
culations of the (post-reionization) kSZ power spectrum. 
In this work, we measure the shape and evolution of 
W 2 (k) in hydrodynamic simulations and apply the re- 
sults to our model for the power spectrum. 

2.4. Final Model 

Our full model for the kSZ power spectrum is given 
by Equation [8l where the expression for including 
non-linear corrections, the gas window function and halo 
contraction is 

^ = ^a 2 f J ||^ 2 (|k - k'|)$S(| k - k'|) 

xP? s L (\k-k'\)p ss (k')i(k,k') ■ (20) 

We remind the reader that the superscript NL represents 
the non-linear matter powe r spectrum, calcula ted using 
the HaloFit procedure of iSmith et al.l (|2003f ). Other- 
wise, we use the linear-theory matter power spectrum. 

One of the principal aims of this work is to explore the 
impact of astrophysical processes such as radiative cool- 
ing, star-formation and supernova feedback on the kSZ 
power spectrum. This is achieved by measuring the win- 
dow function W 2 (k) in hydrodynamic simulations both 
including and omitting these processes. By plugging the 
results into our model, we are able to calculate the effect 
of these processes on the kSZ power spectrum. To inves- 
tigate the impact of non-linear structure formation and 
halo contraction, we measure the curl component of the 
momentum power spectrum, , directly from our sim- 
ulations and compare with the analytic prediction given 
by Equation [20j 



Based on our hydrodynamic simulations, we investi- 
gate three models for the kinetic SZ power spectrum, 
labeled DM (dark matter), NR (non-radiative) and CSF 
(cooling and star- formation) . The three models differ 
principally in the gas window function, W 2 {k), used to 
relate the non-linear dark matter power spectrum to the 
gas density power spectrum, where i = {1, W^ R , W^ SF } 
for the {DM, NR, CSF} models. For the CSF model we 
also include the halo contraction correction to the dark 
matter density power spectrum, $(fc,z), given by Equa- 
tion [TU For the DM and NR models we set this equal to 
1 at all fc and z. 

In our model, simulations are used to determine the 
gas density power spectrum whereas the velocity modes 
are entirely calculated from linear theory. Our technique 
thus circumvents the problems relating to the truncation 
of large-scale velocity modes due to a limited simulation 
box size (see Section W% . Indeed, our simulations accu- 
rately resolve the range of spatial scales over which den- 
sity fluctuations contribute significantly to the kinetic SZ 
power spectrum for 1000 < £ < 10000. We now demon- 
strate this point explicitly. 

In FigureQ]we plot the contribution to the total kinetic 
SZ power at £ = 3000 by density fluctuations at a given 
rcdshift and spatial scale k = 3000/x. In each panel, the 
contours enclose the regions that contribute 25, 50 and 
75% of the total power (inner to outer contours). The top 
panel shows contours for the full kinetic SZ power spec- 
trum including the non-linear density fluctuations. The 
lower panel shows the same contours for the Ostriker- 
Vishniac effect, in which only the linear- regime matter 
power spectrum is used. The dashed blue lines show the 
relation k = £/x(z) at £ = 3000. The shape of the con- 
tours follow that expected for the distance-size relation 
at a constant angular scale. In this figure W 2 (k) is fixed 
to 1 at all scales. 

Comparing the two panels clearly demonstrates the 
impact of non-linear density perturbations on the kSZ 
power. The full kSZ contours extend to larger k and 
lower redshift, where the impact of the non-linear cor- 
rections to the power spectrum are largest. More than 
50% of the full kSZ power is sourced by density fluctua- 
tions in the range 0.6 < k < 5 /iMpc -1 and < z < 3, 
while wavenumbers up to k = 8 /iMpc -1 provide a 
non- negligible contribution. In contrast, we find that 
the principal contribution of velocity modes to the kSZ 
signal comes from the range of scales encompassed by 
0.005 < k < 0.5 /iMpc -1 and thus from much smaller 
k than the density fluctuations. This emphasizes that 
the kinetic SZ power spectrum is generated by small 
scale density fluctuations caught up in large scale veloc- 
ity flows. The OV signal is spread over a larger redshift 
range, but limited to contributions from density fluctu- 
ations at fc < 2 /iMpc -1 . Our simulations accurately 
resolve density fluctuations over the full range of fc plot- 
ted and are thus well-suited to calibrating our model. 

3. SIMULATIONS 

Our simulations are performed using the Adaptive 
Refinement Tree (ART) TV-body+ga s-dynamics code 
(jKravtsov et alJl2002b iRudd et aljl2008[ ). which is an Eu- 
lerian code that uses adaptive refinement in space and 
time, and non-adaptive refinement in mass ([Klypin et al.l 
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TABLE l 
List of Simulations 



Name 


^box 

(Ti-iMpc) 




n A 


li 


Q b h 2 


<?8 


n s 




m v 

(10 9 ft- 1 M Q ) 


Ax 
(h-ikpc) 


L60NR 


60 


0.3 


0.7 


0.7 


0.021 


0.9 


1.0 


256 3 


0.92 


1.8 


L60CSF 


60 


0.3 


0.7 


0.7 


0.021 


0.9 


1.0 


256 3 


0.92 


3.6 


BolshoiNR 


250 


0.27 


0.73 


0.7 


0.023 


0.82 


0.95 


1024 3 


1.08 


3.8 



Note. — ^box is the simulation box side-length, noM and m p are the number of dark of dark matter particles and their mass (in units 
of 10 9 /i _1 Af0), and Ax is the peak spatial resolution of the simulation. 



2001) to achieve the dynamic range necessary to resolve 
the cores of halos formed in self-consistent cosmological 
simulations. 

The simulation properties, their associated box sizes, 
cosmological parameters, and resolutions are summa- 
rized in Table [T] Two of the three simulations (Bol- 
shoiNR and L60NR) are performed in the non-radiative 
regime where baryons are shock heated during structure 
formation but unable to radiatively cool. BolshoiNR is 
a large simulation designed to study the formation of 
galaxy clusters. L60NR and L60CSF are realizations of 
the same set of initial density fluctuations simulated with 
varying physical processes included and were previously 
used to study baryon effec ts on the matter power spec- 
trum in IRudd et al.l (j2008l) . 

A complete description of the physical processes im- 
plemented in the L60CSF simulation can be found in 
IRudd et al.l ([2008). Briefly, gas is converted to colli- 
sionless stellar particles using a Kennicutt- Schmidt type 
density relation for a gas consumption timescale of ~ 
2.7Gyr. Mass, energy, and metals are returned to the 
ISM through prescriptions for types II and la supernovae 
and stellar mass-loss via winds. Radiative cooling rates 
are tabulated over a range of gas temperature, density, 
and mctallicity includ ing a redshift-dependent cosmolog- 
ical UV background (IHaardt fc Madaui ri996) using the 
CLOUDY code (ver. 96b4; iFerland et al.lTl99a) . In con- 
structing these rates the collisional and UV ionization 
equilibria are explicitly calculated. This allows the di- 
rect determination of the fraction of ionized electrons in 
each mesh cell, x, rather than assuming a universally 
averaged value. 

One of the main effects of cooling and star forma- 
tion is to lower the volume-averaged gas density sig- 
nificantly below the cosmic baryon density due to the 
conversion of gas into stars. At z — 4 the gas fraction 
is Pg/pbai = 0.99 (where pbar is the mean baryon den- 
sity of the Universe), whereas by z = it has fallen 
to 0.8. As noted by IRudd et al.l (|2008l) . the conver- 
sion of gas into stars in the simulation is significantly 
more efficient than is observed - especially in the high- 
est density region s - due to the well-kn o wn ov er-coolinq 
problem (see also iKravtsov et al.l [20051 l2009f ). At the 
high-mass end, the halos in this simulation have a stel- 
lar mass fraction at z — that is approximately 50% 
great er than observed for group and cluster -mass ob- 
jects (|Gonzalez et al.ll2007tlGiodini et al.ll2009fl . We note 
that, while the star-formation rate at high redshift z > 4 
is likely to be underestima ted (as many halos are poorl y 
resolved at these epochs, ISpringel &; Hernquistl 120031 ). 
Figure [T] indicates that the majority of the kSZ power 
at £ = 3000 is sourced at lower redshift. Overall, the 



L60CSF simulations should therefore underestimate the 
amplitude of the gas density - and thus the kinetic SZ - 
power spectrum. On the other hand, the absence of these 
processes in our non-radiative simulations results in a 
mean gas density in halos that is substantially larger than 
that observed. The BolshoiNR simulation should thus 
over-estimate the kinetic SZ power. Therefore, while nei- 
ther simulation represents the real Universe, we expect 
them to provide well-motivated lower and upper limits to 
theoretical estimates of the kSZ power spectrum. These 
limits can be utilized to place constraints on the kSZ 
signal from reionization. 

To compare with these two limiting cases we also reran 
the L60CSF simulation having turned off radiative cool- 
ing at z = 1.8 (henceforth referred to as L60CSFz2). 
While this run should be considered unrealistic, it serves 
two useful purposes. Firstly, it allows us to investigate 
whether a radical change in the star-formation prescrip- 
tion can significantly modify the shape of the window 
function, W 2 (k), and the kSZ power spectrum away from 
that predicted by either the BolshoiNR or L60CSF simu- 
lations. This then enables us to evaluate how robust our 
estimate of W 2 (k) is to changes in baryon physics: does 
it vary smoothly between the non-radiative and over- 
cooling case, or does it depend sensitively on the star- 
formation prescription? Secondly, the final stellar (and 
gas) mass fraction in halos of mass M500 > 10 12 h~ x M© 
at z = in this simulation lies between the values pro- 
duced by the BolshoiNR and L60CSF simulations. For 
example, for halos of mass M500 = 5 x 10 13 M Q , the ra- 
tio of stellar (gas) mass to total mass in the L60CSFz2 
simulation is approximately 0.12 (0.8) of the universal 
baryon fraction, compared to (0.9) for BolshoiNR and 
0.6 (0.45) for L60CSF. Therefore, this run provides a 
useful intermediary case between the full cooling plus 
star- formation and non-radiative cases. In Section [5] we 
discuss future work that is required to make precision 
estimates of the kinetic SZ power spectrum. 

The BolshoiNR and L60CSF are the primary simula- 
tions used to gauge the effects of cooling and star forma- 
tion and to calibrate our models. As indicated in Tabled] 
the BolshoiNR simulation has the largest box size (250 
Mpc//i). The range of scales resolved by this simulation 
encompasses the wavenumbers 0.03 < k < 100 ft,Mpc _1 
(we conservatively use an upper limit of 1 /8 the Nyquist 
wavenumber corresponding to the spatial resolution). 
Mildly overdense structures (S « 10) are followed at 
a grid resolution of approximately 60 /ikpc" , corre- 
sponding to an upper limit on k of about 7 /iMpc -1 . 
Therefore, this simulation has adequate resolution in 
low density regions to account for the potential con- 
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tribution of filamentary str uctures to the kSZ signal 
(|Atrio-Barandela et al.l 12008ft . The L60CSF simulation 
has a box size of 60 Mpc/h, resolving fluctuations over 
the range 0.1 < fc < 110 /iMpc" 1 . The resolution in 
lower-density regions is approximately 200 /ikpc -1 , so it 
is possible that density fluctuatio ns in these regions are 
slightly underestimated. However, [Hallman ct al.l (|2009T > 
demonstrate that the contribution of regions of density 
S < 50 to the kSZ power spectrum in their simulations 
is an order of magnitude below that of denser regions. 
Therefore, comparing with the upper panel of Figure [I] 
it is clear that our simulations adequately resolve density 
fluctuations over the range of scales that contribute the 
bulk of the kSZ signal. 

While we do not directly present results from the 
L60NR simulation, we use this to check the effect of vary- 
ing the simulation box size (by comparing BolshoiNR and 
L60NR) and separate it from the effect of adding bary- 
onic physics (by comparing L60NR and L60CSF). 

4. RESULTS 

4.1. Gas Window Function 

In the left panel of Figure [2] we plot the window func- 
tion W 2 (k) = P g as(fc)/ft)M(fc) for the BolshoiNR simu- 
lation, where P gas and Pdm are the measured gas and 
dark matter density power spectra. The points show 
the results over a range of timesteps, corresponding to 
a = 0.4 to 1.0 in steps of 0.1 (from upper to lower at 

fc = 2 /iMpc -1 ). The black solid lines are fits to the 
simulation results using the fitting function given below. 
The red lines show th e linear-theory window function of 
IGnedin fc Hull <|1998f ) given by Equation IT71 

Fluctuations in the gas density follow those of the dark 
matter at large scales (fc < 1 ft-Mpc -1 ), but are rapidly 
suppressed towards smaller scales as the gas thermal 
pressure begins to counter gravity in overdense regions. 
Note that the largest value of k plotted is significantly 
below the maximum that is resolvable for the simulation. 
Hence, the truncation of power at small scales in P gas is 
due to a physical smoothing of the gas density by ther- 
mal interactions rather than an artificial smoothing due 
to the finite resolution of the simulation. 

We find that the following fitting-function provides a 
good description of the BolshoiNR simulation results, 

a) = 0.5 (exp (-k/k' f ) + x + j fc/fc , )7/2 ) , 

(21) 

where k' f = 12.6a" 1 + 6.3 and g(a) = 0.84a" 1 . We have 
verified that Equation [21] also provides a good match to 
the window function measured from the L60NR simula- 
tion. The difference between the BolshoiNR and L60NR 
simulation box sizes and cosmological parameters do not 
significantly affect our results. 

At the final output, the characteristic filter scale, k'j, is 
approximately a factor of 3.3 less than that of the linear- 
theory prediction given by Equation [18] This is also 
evi dent from com p aring the window function suggested 
by IGnedin fc Huil (|1998l red lines) with the simulation 
results; the smoothing of the gas density fluctuations 
occurs at much smaller sca les than in the simul ation. 
The difference between the IGnedin fc Huil ([1998) win- 



dow function and the simulation results is due to shock 
heating of the gas in the simulation as density pertur- 
bations become non-linear. To calculate the Jean's scale 
in Equation [18] we assume that the gas is initially cou- 
pled to the CMB temperature, but then evolves adia- 
batically, T gas oc 1/a 2 , until it reaches a minimum tem- 
perature of 300K (which is also imposed in the simula- 
tion). The lower gas temperature at late times results 
in a large Jeans wavenumber, and thus the smoothing 
of the gas density is limited to very small scales (e.g., 
{kj,kf} w {32,63} /iMpc" 1 at a = 1). However, in the 
simulation, as halos collapse and grow, accretion shocks 
heat much of the gas to temperatures significantly above 
the 300K temperature floor (note that this simulation, 
unlike the L60CSF run, does not include a UV back- 
ground). The mean (mass- weighted) gas temperature in 
the simulation is in excess of 10 6 K at z = 0. Thus we see 
smoothing of the gas distribution at much larger scales 
in the simulation than predicted by purely adiabatic evo- 
lution. 

On the right side of Figure [5] we show the window func- 
tion, WcsF(k), measured in our L60CSF simulation at 
the timesteps corresponding to a = 0.3 to 1.0 in steps of 
0.1 (black lines, from top to bottom at k = 1 /iMpc -1 ; 
for clarity we represent the simulation results with lines 
rather than symbols). In this simulation we are able to 
measure the ionization fraction of hydrogen in each cell. 
We therefore plot the ratio of the free-electron (rather 
than the gas) density power spectrum to that of the dark 
matter. However, since the global neutral fraction is very 
small, using the gas density produces similar results. 

The inclusion of gas cooling and star-formation pro- 
duces a very different window function to the non- 
radiative simulation. Firstly, the overall amplitude 
at large scales is clearly time-dependent. At k = 
0.2 /iMpc" 1 , W 2 {k) varies from 0.92 at a = 0.2 to 0.62 
at a = 1.0. For the BolshoiNR simulation we found that 
W^ R (k — 0.2) re 1 at every timestep. The difference is 
due to the evolving gas fraction within halos in the CSF 
simulation. At a = 1, the square of the ratio of the gas 
to total baryon density in the box is (p g //?bar) 2 = 0.64 
(shown as the horizontal dot-dashed line). The low gas 
density within halos in this simulation significantly re- 
duces the power spectrum of gas density fluctuations 
relative to the dark matter. This is particularly evi- 
dent at low k as the two-halo contribution to the gas 
density power spectrum is prop ortional to the square o f 
the mean gas density in halos (Sembolo ni et al.ll2011|) . 
Therefore, as the gas density is reduced, so is the ampli- 
tude of W^ SF (k) at large scales. 

At intermediate scales W 2 (k) slowly decreases, but 
starts rising again at fc > 5 /iMpc" 1 . The increase in 
power towards small scales is due to the cooling of gas 
at the center of halos, steepening the density profile at 
small radii relative to the dark matter and thus boosting 
the small-scale power. For earlier outputs (upper lines) 
we see that the power in the free-electron density can 
greatly exceed that of the dark matter (Wcgp(fe) > 1). 
With time, this bias gradually disappears due to two ef- 
fects. Firstly, as the simulation progresses, the cold gas 
in halo cores is converted to stars. Secondly, the gas in 
the very center of halos can cool to the point at which 
hydrogen is no longer ionized, reducing the free-electron 
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Fig. 2. — (left) The ratio of the gas density to dark matt er density powe r spect ra, W 2 (k), measured in the BolshoiNR simulation (black 
symbols) and the linear perturbation theory prediction of Gncdin & Hui (1998) (red solid lines). For each, the results are shown from 
a = 0.4 to 1 (where a is the scale factor for the size of the universe) in step s of Aa = 0.1 (from top to bottom at k = 2 hMpc" 1 ). The 
black solid lines show the fit to the simulation points given by Equation 1211 (right) The ratio of the free electron density to dark matter 
density power spectra in the L60CSF simulation (black lines). The results are shown from a = 0.3 to 1 in steps of Aa = 0.1 (from top 
to bottom). Also plotted is the ratio of the free electron density to dark matter density power spectrum for the L60CSFz2 simulation (in 
which radiative cooling is turned off at z = 1.8) at a = 0.4, 0.5, 0.68 and 1 (from top to bottom at k = 10 h Mpc -1 ) 



density. However, the neutral fraction always remains 
small and thus the former is the dominant effect. Us- 
ing the gas (rather than free-electron) power spectrum 
produces very similar results. 

The blue lines in the right-hand panel of Figure [2] show 
the window function obtained for the L60CSFz2 simula- 
tion, in which radiative cooling is artificially turned off at 
z s — 1.8. The results are given for the timesteps corre- 
sponding to a = 0.4,0.5,0.68 and 1 (from top to bottom 
at k = 10 /iMpc -1 ). Switching off radiative cooling has 
a clear effect on the gas distribution across the range of 
scales probed. After z a s , the cold, dense clumps of gas in 
halo cores begin to heat and expand, rapidly suppressing 
small-scale density fluctuations. The reduction of small- 
scale power thus occurs in both CSF simulations, but for 
very different reasons: in L60CSF cold gas is converted 
into stars, whereas in L60CSFz2 the gas mixes with its 
warmer surroundings and begins to expand. This also 
results in a different behavior at large spatial scales. In 
L60CSFz2, the gas becomes less centrally concentrated 
in halos, reducing the suppression in power (relative to 
the dark matter) that is seen in the L60CSF simulation. 
As the simulation progresses towards z = 0, W 2 (k) thus 
increases slightly at small k. By the final output, the 
shape of the window function is beginning to resemble 
that of the BolshoiNR simulation (left-hand panel of Fig- 
ure [2]). Therefore, once cooling is turned off, the window 
function obtained from the L60CSFz2 simulation begins 
to transition from a CSF-like shape to a NR-like shape. 
We show the effect of this on the kSZ power spectrum in 
Section [5] 

The form of W / csf(^ c ) * s evidently more complicated 
than W^ R (fc). To incorporate the CSF window functions 
into our calculation of the kinetic SZ power spectrum, 
we instead interpolate between the simulation results for 
the redshifts between outputs. We assume that W 2 (k) 



remains fixed at its value at k = 0.2 /iMpc~ for scales 
with k < 0.2 /iMpc -1 . For epochs earlier than a = 0.3 
(z > 2.33), we assume that WQ SF (k) smoothly converges 
towards unity at all scales. 

4.2. Momentum Power Spectrum 

We now compare the power spectrum of the curl com- 
ponent of the gas momentum field, A^, predicted by our 
model with that measured directly from the simulations. 
This comparison enables us to test whether our model 
accurately captures the impact of non-linear structure 
growth and baryonic physics in the simulations. This is 
of particularly importance as the kSZ power spectrum is 
a weighted integral of over redshift. We also use our 
model to demonstrate the impact on A^ of the trunca- 
tion of large-scale velocity modes due to the finite size of 
our simulated volumes. 

To measure A^ in our simulations we use a 1024 3 mesh 
rather than the fully refined mesh in each simulation. 
Hence, we are limited to a narrower range of scales than 
for the density power spectra. For the BolshoiNR simu- 
lation (which has a box side-length of 250 Mpc/h) this 
range is 0.03 < k < 12 ftMpc" 1 . For the L60CSF simu- 
lation (60 Mpc/h box) the range is 0.1 < k < 50. 

In Figure [3] we plot A^ for the BolshoiNR simulation 
at four timesteps, corresponding to z = 4, 1.87, 1 and 0. 
The momentum power is shown in terms of the dimen- 
sionless quantity Ask/ 'H(z). The black dots represent 
the measured simulation power spectrum. The blue lines 
give the analytic prediction in two regimes: the linear 
regime calculation (the Ostriker-Vishniac effect) given 
by Equation [TT] (dotted), and our fiducial model given 
in Equation l201 which includes the non-linear corrections 
to the density power spectrum (solid). Both include the 
window function, W^ R (k). 

Non-linear density fluctuations become significant 
when Af,k/H(z) fts 1. Our full non-linear model re- 
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Fig. 3. — The power spectrum of the curl component of the momentum field in the BolshoiNR simulation at four timesteps: z = 4, 
1.87, 1 and (top left, top right, lower left, lower right). The black dots represent the si mu lation results, the blue lines give the analytic 
predictio n of Ask/ H(z) in two regimes: the linear regime calculation given by Equation [Tl] (dotted) , and the non-linear version given by 
Equation [13] (solid). The red line shows the 'truncated' model (see text). 



produces the simulation results extremely well at k > 
0.1 /iMpc -1 , but lies systematically above the simula- 
tion points at lower k. This is simply due to the finite 
volume of the simulation box. A 2 B (k) is generated by a 
convolution between velocity modes at wavenumber k! 
with density modes at \k — k\. The largest contribution 
comes from velocity modes at large scales coupling to 
density modes at smaller scales. Velocity modes on scales 
larger than the simulation box size are not accounted for, 
resulting in an underestimate of the momentum power. 

To demonstrate that this is indeed the case, we re- 
calculate Ag, placing an upper and lower limit on the 
integral over k 1 in Equation [201 such that only the ve- 
locity modes encompassed by the simulation box are in- 
cluded. The red lines in each panel of Figure [3] show the 
results of this 'truncated' model. It accurately repro- 
duces the simula tion on all sc a les. T his was previously 
demonstrated by iZhang et al.l (|2004j ). who performed a 
similar test by calculating a discretized version of Equa- 
tion [20] and compared the results with their simulations. 
Despite the overall very good agreement, the truncated 
model does slightly underestimate the simulation results 
at k > 2 /iMpc^ 1 . This may be due to the impact of 
a non-linear, curl component, of the velocity field gen- 
erated by shell-crossing dur i ng halo formatio n (as sug- 
gested by [Zhang ct al. 2004; Sh ao et all 120111 ). To mit- 
igate this, Z hang et all (|2004D suggest replacing P$s in 



Equation [T3] with it's non-linear counterpart. When we 
do so, we find that the ratio of model to simulation in- 
creases slightly, from 0.88 to 0.94 at k = 2 /iMpc -1 . 

In Figure @] we plot the curl component of the gas 
momentum power spectrum in the L60CSF simulation 
(black dots) at the same four timesteps. In this figure, 
the blue lines show the results omitting (dashed) and in- 
cluding (solid) the halo contraction correction applied 
to the (non-linear) density power spectrum described 
in Section 12.21 Both lines include the window function 
WQ SF (k) measured in the L60CSF simulation. 

At the low-fc end, the simulation points again lie sys- 
tematically below the model. However, for this simu- 
lation the discrepancy extends to much higher k. The 
L60CSF simulation box is more than a factor of four 
times smaller than the BolshoiNR simulation (see Ta- 
ble [1]), and so the velocity mode truncation is more se- 
vere. We note that our L60NR simulation underesti- 
mates (relative to our non-radiative model) by a sim- 
ilar amount, confirming that the box-size effect is inde- 
pendent of the baryonic physics included in the simula- 
tion. As in Figure[3J the red lines show our model predic- 
tion having limited the calculation to include only those 
scales encompassed by the simulation. The model again 
reproduces the simulation results very well (although at 
a = 1 it appears to slightly overestimate the momentum 
power at scales around k = 1 /iMpc -1 ). 
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Fig. 4. — Similar to Figure[3]but for the L60CSF simulation. The black dots represent the simulation results, the solid and dashed blue 
lines give the analytic prediction of Agk/H(z) including and omitting the halo contraction correction, respectively. The red line shows the 
'truncated' model (see text). 



The dashed lines demonstrate the model with no halo 
contraction correction. It is clear that, without this cor- 
rection, the model underestimates for wavenumbers 
k > 10 hMpc~~ , especially at z = 0. The solid blue line 
demonstrates that our simple correction does a reason- 
able job of reproducing the simulated momentum power 
spectrum. However, as Figure [T] demonstrates, these 
scales do not contribute significantly to the kSZ power 
spectrum at I = 3000. 

5. KINETIC SZ POWER SPECTRUM 

We now utilize the gas window functions measured in 
our NR and CSF simulations to explore the impact of 
baryonic physics on the shape and amplitude of the ki- 
netic SZ power spectrum. As described in Section 12.41 
we evaluate three baseline models: DM, NR and CSF. 
For the DM matter model W 2 (k) = 1, for the NR and 
CSF we use the window functions W^ R csfW measured 
in our simulations. We also investigate the redshift de- 
pendence of the kSZ power spectrum, and its scaling 
with cosmological parameters for each model. Finally 
we compare the results presented in this work with those 
of previous theoretical studies as well as the latest obser- 
vational constraints. 

5.1. The Impact of Baryon Physics 

In Figure [S] we plot the kinetic SZ power spectrum for 
our three models: DM, NR and CSF (solid black, red 



and blue lines). The dashed lines represent variants to 
the DM and CSF models. The black dashed line shows 
the linear theory version of the DM model, which we la- 
bel OV (Ostriker-Vishniac). The blue dashed line shows 
the kSZ power spectrum calculated using the L60CSFz2 
simulation window function. The power spectra are plot- 
ted in terms of D e = 1(1 + 1)Q/2tt. 

Comparing the dashed and solid black lines demon- 
strates the impact of non-linear structure growth on the 
kinetic SZ power spectrum. Including the non-linear cor- 
rections significantly boosts kSZ power at all but the 
largest angular scales. At £ = 3000, non-linear correc- 
tions have increased the kSZ signal by a factor of 2. By 
£ = 10, 000, this has increased to a factor of 3.5. The 
kSZ effect is sourced by small-scale density fluctuations 
caught up in large-scale bulk velocity flows. Therefore, 
the large boost in the amplitude of small-scale density 
perturbations due to non-linear gravitational collapse 
also enhances the kSZ power, especially at small angular 
scales. 

The difference between the DM and NR models are 
small. At £ = 3000 (10, 000) the BolshoiNR window func- 
tion reduces the kSZ power from the DM model by 0.20 
(0.72)/iK 2 . The NR window function suppresses power 
at scales k > 8 /iMpc -1 (see the left panel of Figure [2]). 
As demonstrated by Figure [TJ gas density fluctuations at 
these scales contribute little to the kSZ power spectrum 
(at £ = 3000). The application of (k) thus translates 
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into only a small reduction of power (relative to the DM 
model) at the angular scales most sensitively probed by 
current small-scale CMB experiments. 

The window function measured in the L60CSF simu- 
lation has a greater impact on the kSZ power spectrum. 
The large-scale suppression of the gas density power spec- 
trum seen in Figure [2] produces a significant reduction in 
kSZ power at all angular scales. At I — 3000, the power 
is reduced by 1.05^K 2 (32%) and by 1.43^K 2 (36%) at 
I = 10, 000, relative to the NR model. As discussed in 
the previous section, this reduction in power is primarily 
driven by the decreased gas density in halos due to over- 
cooling and excessive star-fo rmation. A s imila r reduction 
in kSZ power was noted by iTrac et ail f|2011h when the 
stellar mass fraction in groups and clusters was increased 
in their semi-analytic model, reducing the gas fraction by 
an equivalent amount. 

The blue dashed shows the results using the window 
function measured in the L60CSFz2 model. The ampli- 
tude of the power spectrum in this case is close to the 
mean of the NR and CSF models, with the shape un- 
changed. A comparison to Figures Q] and [2] provides a 
clear explanation. Figure [1] demonstrates that the peak 
of the kSZ contribution comes from spatial scales corre- 
sponding to 1 /iMpc -1 < k < 3 /iMpc -1 at redshift 
z ~ 1. At these scales and redshifts, the L60CSFz2 
window function lies almost directly between that of the 
BolshoiNR and L60CSF simulations, and thus the kSZ 
power spectrum also lies in between these two cases. The 
dominant factor determining the amplitude of kSZ power 
thus appears to be the mean gas density in halos, which 
is controlled by the star-formation rate. 

We have also investigated the impact of non-linear ve- 
locity fluctuations by replacing the linear density power 
spectrum (i.e. the velocity component of t he A^) with 
its no n-linear counterpart, as suggested by iZhang et al.l 
(2004). We find that, at £ = 3000 (10,000), the ampli- 
tude of the kSZ power spectrum in the NR model in- 
creases by 1% (4%). We have chosen to omit this correc- 
tion from our fiducial models, as the velocity contribu- 
tion to A 2 B in Equation [T3l is determined by the gradient 
of the density field and is therefore curl-free. Hence, it 
is not clear how one should correctly account for non- 
linearities in the velocity field. Nevertheless, the impact 
on the kSZ power spectrum is likely to be small. 

In Figure [5] we plot the relative contribution to the 
kSZ power spectrum of slices in redshift over the range 
< z < z re i, at £ = 3000. The solid lines show the 
differential contribution, dDi/dz, the dashed lines show 
the cumulative contribution, Di(< z)/Dg(< z re i), where 
Zrei — 10. The line colors represent our different mod- 
els, DM (black), NR (red), CSF (blue) and also the OV 
contribution to the DM model (grey). For comparison, 
we also plot the results obtained using the L60CSFz2 
window function (green). 

Comparing the differential redshift contribution of the 
DM and OV models again demonstrates the enhance- 
ment provided by non-linear structure formation at low 
redshift. Half of the kSZ signal in the DM model comes 
from z < 2, whereas the equivalent fraction for the OV 
is attained at much higher redshifts (z ~ 5). 

The NR model has a similar redshift distribution to 
the DM model. However, the CSF model predicts a red- 
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Fig. 5. — The kinetic SZ power spectrum for our three models: 
DM, NR and CSF (solid black, red and blue lines). The black 
dashed line shows the linear theory version of the DM model (the 
Ostriker-Vishniac effect). The blue dashed line shows kSZ power 
spectrum calculated using the L60CSFz2 window function. The 
vertical dotted lines denote the approximate range in I over which 
current small-scale CMB experiments such as SPT and ACT are 
sensitive to kSZ power. 

shift distribution that is significantly flatter. While the 
differential redshift distribution peaks at the same red- 
shift, the amplitude of the peak is considerably lower. 
The half-way point of the cumulative distribution is at 
a much higher redshift, z ~ 4, than for the DM and 
NR models. The right panel of Figure [5] shows that the 
effects of cooling and star-formation, which reduce the 
gas density in halos, become increasingly more signifi- 
cant towards lower redshift. This counteracts the boost 
to the kSZ signal provided by non-linear density fluctua- 
tions such that the redshift distribution of the CSF more 
closely resembles that of the OV model. 

The L60CSFz2 simulation, in which radiative cooling 
is turned-off at z = 1.8, produces a kSZ redshift distri- 
bution that lies almost exactly in between the NR and 
CSF models. The distribution diverges from the L60CSF 
line at z = 1.8 and then remains at roughly 2/3 the am- 
plitude of the NR model. This is to be expected; once 
cooling is turned off, gas is no longer being consumed 
by star-formation. Hence, at lower redshifts, the differ- 
ences between the L60CSFz2 and the NR model cease to 
grow. These results suggest that other (and less sudden) 
variations in the star-formation history would produce a 
predictable change in the kSZ redshift distribution: the 
decrease compared to the NR case would depend on the 
integrated star-formation rate to that redshift. While 
it may underestimate the star-formation rate at high 
redshift (due to the limited mass resolution), the CSF 
model only diverges significantly from the OV prediction 
at z < 4. Given the over-cooling problem at low red- 
shift, it is reasonable to assume that the NR and CSF 
cases encompass the expected redshift distribution of the 
kSZ power spectrum at £ = 3000. 

The redshift distribution of the kSZ power spectrum 
can potentially be measured usi ng kSZ tomography 
(|Ho et al.l [20091: iShao et al.l [20Tlh . The tomography 
method utilizes catalogues of galaxies with precisely mea- 
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Fig. 6. — The contribution of slices in redshift to the kSZ power 
spectrum at i = 3000. The solid lines show the differential con- 
tribution, the dashed lines show the cumulative contribution. The 
line colours represent our different models, DM (black), NR (red), 
CSF (blue), and the OV contribution to the DM model (grey). Also 
plotted is the kSZ redshift distribution predicted by the L60CSFz2 
simulation (green). 

sured spectroscopic redshifts to reconstruct the large- 
scale velocity field from the 3-D galaxy distribution. 
When integrated along the line of sight, the product of 
the reconstructed velocity and density fields provides an 
estimator for the kinetic SZ temperature fluctuations. 
This estimator can be broken down into redshift slices 
and cross-correlated with a map at CMB frequencies en- 
abling the true kSZ contribution from within that red- 
s hift slice to b e mea sured. 

iShao et ail ([201 1[ ) de monstrate that, by combining 
Planck with BigBOSS (|Schlegel et all 120111 ) it will be 
possible to measure the kSZ power from the redshift 
range 0.2 < z < 0.6 to better than 10% precision at 
£ = 3000, and within the range 0.6 < z < 1 to better 
than 20%. Utilizing data from higher angular resolution 
CMB experiments such as ACT should significantly in- 
crease the precision of these measurements. This will not 
only provide a unique, temperature- independent probe of 
the inter-galactic medium at lower densities than those 
typically probed by X-ray and tSZ surveys, but would 
place tight constraints on the models of the kSZ power 
spectrum presented here. Figure [6] demonstrates that 
much of the theoretical uncertainty on the kSZ signal 
comes from within the range < z < 2; kSZ tomogra- 
phy can provide a powerful probe of the kinetic SZ effect 
over the lower half of this redshift range. 

5.2. Cosmological scaling 

We have discussed in detail the theoretical uncertainty 
in the kinetic SZ power spectrum from astrophysical 
sources. However, the kSZ signal also scales sensitively 
with cosmological parameters. There is therefore an ad- 
ditional uncertainty on theoretical estimates of the kSZ 
power spectrum due to the precision with which these 
parameters have been measured. In this section we in- 
vestigate how the amplitude of the kSZ power spectrum 
scales with cosmological parameters - in particular, erg 



and t - and how this scaling varies between our NR and 
CSF models. 

We assume that scales as a power-law with each 
cosmological parameter p, i.e., Di oc p a . We vary each 
parameter by 20% of its fiducial value (holding the other 
parameters fixed) and measure the value of a. We have 
verified that the power-law scaling is appropriate for all 
parameters. The results are presented in Tables [2] and [3J 
We give the results over the angular multipole number 
range 1000 < £ < 10,000. Also given are the band- 
powers, De, for each model. 

The kSZ power spectrum amplitude is most sensitive 
to cr 8 , scaling as of a ~ 51 from £ = 1000 to 10,000 (NR 
model). This scaling is simple to understand. The kSZ 
power spectrum involves the product of the density and 
velocity power spectra; in the linear regime each con- 
tributes a power of 2 to the scaling with erg. The OV 
power spectrum thus scales exactly as cr|. Non-linear 
corrections to the matter power spectrum steepen the erg 
dependence. This steepening is related to the fraction of 
the kSZ signal that originates at low redshift where the 
non-linear corrections have grown large. Figure [6] shows 
that the maximal contribution to the OV signal comes 
from z = 2, whereas the full kSZ signal for the DM model 
peaks at z — 0.5. This shift in the redshift contribution 
is driven by non-linear structure formation, which ex- 
plains the steeper erg dependence. The same argument 
applies to the difference in the erg scaling between the NR 
and CSF models; a smaller fraction of the signal derives 
from low redshift in the latter and so the scaling with erg 
is weaker. We note that changing erg principally varies 
the overall amplitude of the kSZ signal; while the scaling 
does have some ^-dependence, the impact on the shape 
of the power spectrum is small over the parameter space 
explored. 

Some of the tightest constraints on erg to date have 
been obtained from measurements of the primary CMB 
power spectr um and the abundanc e of galaxy clusters. 
For example, iVikhlinin et al.l (|2009f) used a sample of 86 
clusters to obtain cr8(fi M /0.25) - 47 = 0.81 ± 0.01 ± 0.02 
(statistical and systematic errors, respectively). Using 
the CSF model scaling, this translates to an uncertainty 
of approximately ±15% on the kSZ power spectrum at 
£ = 3000. This is equivalent to the uncertainty in Di 
provided by baryon physics presented in the previous sec- 
tion. 

The magnitude of the kSZ signal is also sensitive to the 
optical depth to reionization, r, which is directly related 
to the redshift at which reionization occurs (see Equation 
[2]). In Tables [2] and [3] we give the scaling of T>i with both 
r and z rc i. The CSF model has a greater sensitivity to 
each of these parameters as a slightly larger fraction of 
the signal is generated at high redshift. In principle, 
the post-reionization kSZ power spectrum could be used 
to measure z re j. In practice, the uncertainty on other 
cosmological parameters washes out this information. 

To this point we have ignored the effect of helium reion- 
ization. The amplitude of the kSZ signal also scales as 
the square of the ionization fraction, x 2 , which is depen- 
dent on the ionization state of helium. Our fiducial model 
assumes that helium remains neutral at all epochs. If, in- 
stead, we assume that helium is singly or doubly ionized 
at z ro ; then the amplitude of the kSZ power spectrum 
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will increase by a (^-independent) factor of 1.16 or 1.33, 
due to the increase in the free electron density. A more 
realistic model in which helium i s singly ionized at z = 6 
and doubly ionized at z = 3 fe.g. lFurlanetto &: Ohll2008l) 
would increase the kSZ power by a factor of 1.22 (1.26) 
at £ = 3000 (10,000) relative to our baseline, neutral he- 
lium, model. Hence, the level of uncertainty on the kSZ 
power spectrum due to Helium reionization is equivalent 
to that due to the uncertainty on <jg. 



TABLE 2 

kSZ Cosmological Scaling: NR. 



£ 




H 


0"8 






r 




fiK 2 


70 km/s/Mpc 


0.82 


0.047 


10.0 


0.076 


1000 


1.97 


1.18 


4.31 


2.29 


0.46 


0.32 


2000 


2.82 


1.51 


4.48 


2.18 


0.46 


0.32 


3000 


3.24 


1.68 


4.62 


2.13 


0.43 


0.30 


4000 


3.49 


1.80 


4.73 


2.10 


0.40 


0.28 


5000 


3.66 


1.88 


4.81 


2.08 


0.37 


0.26 


6000 


3.77 


1.95 


4.88 


2.06 


0.34 


0.24 


7000 


3.85 


2.00 


4.94 


2.05 


0.32 


0.22 


8000 


3.90 


2.05 


4.99 


2.03 


0.30 


0.21 


9000 


3.93 


2.09 


5.04 


2.02 


0.29 


0.20 


10000 


3.95 


2.13 


5.07 


2.01 


0.28 


0.19 



Note. - The scaling a of the kSZ power spec- 
trum with cosmological parameters, p, for the NR model, 
where we assume Di oc p a . The fiducial parameters are 
given in the second line and the corresponding values of 
T>i in the second column. 



TABLE 3 
kSZ Cosmological Scaling: CSF 



I 




Ho 


o~8 


fib 


■^rei 


T 






70 km/s/Mpc 


0.82 


0.047 


10.0 


0.076 


1000 


1.43 


1.09 


4.19 


2.31 


0.63 


0.43 


2000 


2.00 


1.46 


4.33 


2.18 


0.66 


0.45 


3000 


2.19 


1.65 


4.46 


2.12 


0.64 


0.44 


4000 


2.27 


1.78 


4.57 


2.09 


0.60 


0.41 


5000 


2.32 


1.87 


4.67 


2.06 


0.55 


0.38 


6000 


2.36 


1.94 


4.76 


2.04 


0.52 


0.35 


7000 


2.40 


2.00 


4.83 


2.02 


0.48 


0.33 


8000 


2.44 


2.06 


4.89 


2.01 


0.45 


0.31 


9000 


2.48 


2.10 


4.95 


2.00 


0.42 


0.29 


10000 


2.52 


2.14 


4.99 


1.99 


0.40 


0.27 



Note. — Similar to Table [2] but for the CSF model. 



5.3. Comparisons with Simulations 

A number of previous studies have used the output 
of cosmological simulations to make predictions for the 
shape and amp litude of the kineti c SZ power spectrum 
(Ida Silva et all 120011: iWhite et all 120021: lHallman et al.l 
I2009t iTrac et al.l 120111: IBattaglia et al.l I2010D . Rather 
than calibrating an analytic model to their simulations 
(as we have done), these studies generate mock sky maps 



of the temperature fluctuations sourced by the kinetic SZ 
effect. Maps are typically constructed by stacking out- 
puts of the simulated volume over a range of timesteps, 
often rotating or translating each output to prevent the 
repetition of structures along the line of sight and gener- 
ate a larger simulated sky area. 

Comparing our model predictions with those made 
from synthetic sky maps is a non-trivial task. The sim- 
ulations encompass a range of cosmological parameters 
and often have a more limited redshift range than our 
fiducial model (where z rG ; = 10; for simplicity we re- 
fer to the highest redshift output of each simulation as 
l z rc i). Furthermore, as we described in detail in Section 
[U the limited size of the simulation volume results in 
the truncation of large-scale velocity modes and thus a 
significant underestimate of the kSZ signal. For exam- 
ple, we estimate that a simulation box of side length 100 
Mpc/h would underestimate the kSZ power at £ = 3000 
by 60%. A box size of at least 1 Gpc/h is required to fully 
account (to within 1% of the total power) for all the ve- 
locity modes that contribute to the kSZ power spectrum 
at £ = 3000. 

To compare our model with that of previous studies, we 
must correct their results to account for the limited sim- 
ulation box size, the variations in the maximum redshift 
assumed, and differing cosmological parameters. We use 
our NR model as a baseline to rescale the simulation 
results to our fiducial cosmology (namely, cr§ = 0.82, 
fib = 0.047 and z rc j = 10), and to estimate the amount 
by which the kSZ amplitude must be increased to ac- 
count for velocity mode-truncation. We also include a 
helium correction such that the level of helium ioniza- 
tion is consistent with our fiducial model (i.e., neutral at 
all epochs). 

The predictions for the kSZ effect from simulations pre- 
sented in previous work are shown in Table HI The ta- 
ble shows three columns for the kSZ power at £ = 3000 

3000 ) j the first shows the raw prediction for D3000 taken 
from each work. The second column shows the kSZ power 
rescaled to our fiducial cosmology using the scalings given 
in Table [2] The third column shows 2?3ooo having addi- 
tionally corrected for the simulation box size in each case. 

We consider the results from the smoothed-particle- 
hydrodynamics (SPH) simulations of W hite et al.l (2002, 
henceforth, WHS02) the Eulerian ENZO simulations of 
lHallman et "ail (|2009l H09) and both the non-radiative 
(NR) simulations and those including cooling, star- 
formation an d active galactic nucle us (AGN) feedback 
performed by IBattaglia et all (pOlOl B10). The WHS02 
simulation includes radiative cooling, star-formation and 
galactic winds (although they note these processes may 
be inhibited by the limited mass resolution of their sim- 
ulation) . The H09 simulation is run in the non-radiative 
regime. We also includ e the 'adiabatic' a nd 'standard' 
semi- analytic models of ITrac et al.l (|2011l TBOll) (see 
also ISehgal et al.ll2Q10T) . The adiabatic model assumes 
that gas resides in hydrostatic equilibrium in the poten- 
tial well of dark matter halos identified in their N-body 
simulation. The standard model also includes simple pre- 
scriptions for star-formation and non-gravitational en- 
ergy feedback from supernovae and AGN. 

The cosmological rescaling for the WHS02 and H09 
simulations are fairly large as both these simulations were 
run assuming as = 0.9 and z rc i = 19 and 3, respectively. 



14 



Shaw et al. 



The box-size corrections are largest for WHS02, who used 
a 100 Mpc/h box, and BIO, who generated their kSZ 
predictions using simulation box sizes of 165 Mpc//i. 

In general, the simulations predict between 2.27 and 
3.91 /iK 2 , depending on the level of gas physics included, 
although H09 measure an amplitude of more than twice 
the upper end of this range. They are therefore consis- 
tent with the predictions of our NR (3.2^iK 2 ) and CSF 
(2.2/iK 2 ) models. BIO find that the inclusion of radia- 
tive cooling, star-formation, and AGN feedback reduces 
the kSZ power by 1/iK 2 . The semi-analytic models of 
TBOll lie on the low end of the range of simulations 
discussed here. Their 'adiabatic' model is similar to our 
NR model, while their 'standard' model includes pre- 
scriptions for star-formation and energy feedback. The 
effect of these prescriptions is to reduce the kSZ signal by 
0.4/iK 2 . They also explore the kSZ signal obtained when 
the stellar mass in groups and clusters was increased be- 
yond that of the standard model, finding a further re- 
duction of 0.3/iK 2 . This supports our assertion that the 
CSF model - which predicts the lowest amplitude of all - 
provides a robust lower limit on the kSZ power spectrum 
due to the reduction of the mean gas density in groups 
and clusters via star-formation. 



TABLE 4 
kSZ Predictions from Simulations 



Paper 




©3000 [iiK 2 ] 






Unnorm. 


Cosm Cor. 


Box Cor. 


WHS02 


3.00 


1.57 


3.91 


H09 


7.40 


8.50 


9.45 


BIO, NR 


2.50 


2.42 


4.03 


BIO, AGN 


1.50 


1.45 


2.42 


TBOll, adiabatic 


2.50 


2.70 


2.70 


TBOll, standard 


2.10 


2.27 


2.27 


This work, DM 






3.44 


This work, NR 






3.24 


This work, CSF 






2.19 



Note. - The amplitude of the kSZ power predicted 
by hydrodynamical simulations in previ ous work. W e show 
the results from the SPH simulations of White et al] (120021 , 
WHS02), the ENZO simulations of IHallman et~afl (|2009l . 
H09), both the non-radiative (NR) simulations and those 
including cooling, sta r-formation and AGN feedback from 
iBattaelia et al] (|2010l . B10) and the 'adiabatic' and 'stan- 
d ard' models from the N-body plus semi-analytic approach 
of lTrac et ail (poTll . TBOll). 



5.4. Comparisons with Observations 

To date, the kinetic SZ effect has not yet been de- 
tected observationally. However, both the Atacama Cos- 
mology Telescope and the South Pole Telescope have 
placed upper limits on the amplitude of the kSZ power 
at i = 3000. From 296 deg 2 of data, iDunklev et all 
(|2011f) o btained a 2a upper li mit of D e = 8/iK 2 . More 
recently IShirokoff et al.l (pOllD obtained an upper limit 
of D e = 6.5^K 2 (also 2a) from 210 deg 2 . Therefore, ob- 
servations do not currently constrain any of the models 
that we have presented in this work, although the pre- 



dictions of the simulations of IHallman et al.l (|2009l ) are 
inconsistent with the most recent measurements. 

The main difficulty in measuring kSZ power lies with 
disentangling the signal from both the thermal SZ effect 
and bright foregrounds, primarily dusty, star-forming 
galaxies (DSFG). Unfortunately, the angular shape of 
both the tSZ and DSFG power spectra is similar to that 
of the kSZ signal, so it is difficult to separate these 
signals in ^-space. However, each has a very different 
frequency dependence; multifrequency observations will 
be extremely effective in separating these components. 
Combining SPT or ACT data with Planck and Herschel 
should enable a significant detection of kSZ power in the 
near future. 

6. DISCUSSION & CONCLUSION 

The kinetic SZ (kSZ) power spectrum is generated by 
the coupling between large scale velocity flows and small- 
scale density perturbations. To predict its amplitude and 
shape it is necessary to understand the behavior of the 
power spectrum of gas density fluctuations over the range 
of scales corresponding to0.1</c<10 /iMpc -1 . In this 
work, we have introduced a new model for the kinetic SZ 
power spectrum that accounts for the effect of baryonic 
physics on the power spectrum of gas density fluctuations 
and thus on the kSZ power spectrum. 

To this end, wc defined a window function, W 2 (k) = 
Pgas{k) / PuM(k) , to provide a mapping between dark 
matter and gas density power spectra in our calculations. 
We utilized hydrodynamic simulations - run in both the 
non-radiative regime and including radiative cooling and 
star-formation - to measure the window functions and 
investigated their effect on the kSZ power spectrum. We 
have presented three models for the kSZ power spec- 
trum: DM (dark matter), in which gas density fluctu- 
ations follow those of the dark matter at all scales, NR 
(non-radiative) and CSF (cooling & star-formation) in 
which we use the window functions measured in our hy- 
drodynamic simulations. 

There is only a small difference between the DM and 
NR models. Gas density fluctuations in our non-radiative 
simulations are suppressed at spatial scales smaller than 
those that contribute significantly to the kSZ power spec- 
trum (for I < 10, 000). At I = 3000, the NR model pre- 
dicts V e = 3.24^iK 2 , only 0.20^K 2 below the DM model. 
However, the CSF power spectrum has a significantly 
lower amplitude as well as a flatter shape than cither 
the NR or DM models. At I = 3000, the CSF model 
predicts V e = 2.19/xK 2 , 1.25/xK 2 below the DM model. 
The reduction in power is driven by the decrease in the 
mean gas density in group- and cluster-mass halos due to 
the high-levels of cooling and star-formation in our sim- 
ulation. This in turn reduces the amplitude of the gas 
density power spectrum and thus the kSZ power spec- 
trum. 

To investigate the impact on the shape and amplitude 
of the kSZ power spectrum of variations in star- formation 
history, we rerun our CSF simulation having turned off 
radiative cooling at z = 1.8. We find that the resulting 
window function smoothly evolved from a CSF-like shape 
to a NR-like shape, with the amplitude (at large scales / 
low k) lying midway between the two cases. The result- 
ing kSZ power spectrum was thus also almost directly 
between that of the NR and CSF models. 
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We have argued that our NR and CSF models provide 
reasonable upper and lower limits to the effect of astro- 
physical processes on the kSZ power spectrum. The NR 
model is calibrated from a simulation that included no 
cooling and star-formation and thus no means of reducing 
the gas density. On the other hand, our CSF model is cal- 
ibrated to a simulation that suffers from the over-cooling 
problem, i.e., from excessive cooling and star- formation. 
Hence, we expect our CSF model to underestimate the 
true kSZ signal. Nevertheless, the difference between the 
NR and CSF models at I = 3000 is only 1/xK 2 . Taking 
the mean of these two models, we find that the astro- 
physical uncertainty on the kSZ power spectrum ampli- 
tude is roughly ±20%. This is significantly less than that 
of the thermal SZ power spectrum, for which the cur- 
rent theoretical uncertainty on the amplitude is ±50% 
(IShaw et alJ!2l)Tnt iTrac et alJl20ll . 

There are two caveats to this argument. Firstly, our 
fiducial models assume helium remains neutral at all 
epochs. Singly or doubly ionized helium would increase 
the amplitude of our models by up to 33%, depending on 
the redshifts at which helium reionization occurs. Sec- 
ondly, the amplitude of the kSZ power spectrum is sensi- 
tive to cosmological parameters, namely as, fif>, Hq and 
r (and, equivalently, 2 re ;). We have investigated the scal- 
ing of our models with cosmological parameters, finding 
that £>3ooo oc a^Cl^H^r - 44 for the CSF model. The 
current la uncertainty on a* obta ined from cluster num- 
ber counts (jVikhlinin et al.ll2009f ) translates to a ±15% 
uncertainty on X>3ooo- 

We have compared our models with pred ictions made 
direc tly from hydrody n amic si mulations (| White et aT] 
2002t Mailman et all 120091: iBattaglia et all 120101: 
Trac et al.l 1201 ID . These studies measured the kSZ 
power spectrum using synthetic kSZ maps constructed 
by the stacking of simulation box outputs. A serious 
drawback of this approach is that velocity modes are 
truncated at the scale of the simulation box. For 
boxes of side length <C 1 Gpc/h, this can result in a 
significant underestimate of the kSZ signal. We used our 
model to correct the predictions of the simulations in 
previous studies for this effect, as well as differences in 
cosmological parameters. We then found a reasonably 
good agreement among different simulations results and 
the range predicted by our models. 

In this work we consider only the homogeneous- or 
post-reionization contribution to the kinetic SZ signal. 
It is well known that models of inhomogeneous reion- 
ization, in which different regions of the Universe are 
reionize d at differen t times , pre dict a patchy kin e tic SZ 
signal. iZahn et al.l (|2005ft and iMcQuinn etaU (2005) 
demonstrate that, to first order, the magnitude of the 
patchy signal is dependent on the duration of reioniza- 
tion. Given current technology, the kinetic SZ power 
spectrum thus provides a unique probe of the redshift 
range spanned by the epoch of reionization. However, 
current observations only provide a measure of the sum 
of the patchy and post-reionization kSZ power spectra at 
angular scales around I = 3000. Therefore, in order to 
extract the patchy component, it is important to have a 
good theoretical understanding of the post-reionization 
contribution. In this work we have performed a detailed 
investigation of the theoretical uncertainty on the post- 
reionization kSZ power spectrum due to astrophysical 



processes, cosmological parameters and helium reioniza- 
tion. Adding these in quadrature produces a total uncer- 
tainty of ^30% on the amplitude of the post-reionization 
kSZ power spectrum at £ = 3000. 

Improvements in our understanding of the kSZ power 
spectrum can be made both theoretically and observa- 
tionally. Models and simulations of the kinetic SZ power 
spectrum can be further developed to make improved 
theoretical predictions. The over-cooling problem may 
be mitigated by the inclusion of energy f eedback from 
active galactic nuclei in our s i mulations dSiiacki et al.l 
20071 [Booth fc Schayd [200l IBattaglia et al.l 120101: 
McCarthy et al.l 120101 : iTevssier et al.l 120111 ). Feedback 
from AGN can heat gas sufficiently in high density re- 
gions to slow the loca l cooli ng (and thus star-formation) 
rate. Ida Silva et al.l ()2001l ) demonstrated that the in- 
clusion of a preheating model to their simulations sup- 
pressed star-formation and increased the kSZ power with 
respe ct to simulation s witho ut preheating. On the other 
hand lBattaglia "et al.l ((2010) also demonstrate that AGN 
also heat gas out to large cluster radii, flattening the 
density profile, reducing the gas fraction in groups and 
clusters. The net impact of AGN feedback on the SZ 
power spectrum will depend on the extent to which these 
two effects balance one another. Recent models of AGN 
in simulations have typically been tuned to match local 
(z ~ 0) observations, however, as we have demonstrated, 
a significant contribution of the kSZ comes from larger 
redshifts where the star formation histories of halos are 
relatively unconstrained. Furthermore, higher (sub-kpc) 
resolution simulations may be required to resolve star- 
formation in lower mass halos. 

From an observational perspective, small-scale CMB 
experiments have currently only placed upper limits 
on the amplitude of the kSZ po wer spectrum and do 
not constrain any o f our models ■ Slii mkoir et all 120711 
iDunklev et al.l 120111 ). However, with increasing area 
and frequency coverage, observations will become signifi- 
cantly more constraining in the near future. kSZ tomog- 
raphy has the potential to constrain models of the post- 
reionization signal by measuring its redshift dependence. 
This can be achieved by cross-correlating a kSZ estimator 
constructed from spectrosc o pic galaxy catalo gues with a 
CMB map (|Ho et al.l 12001 IShao et al.l[2lTfl . We have 
investigated the contribution of slices in redshift to the 
total kinetic SZ power for each of our models. We found 
that the principle difference between the NR and CSF 
models comes from sources in the range < z < 2. Up- 
coming surveys have the potential to measure the red- 
shift distribution of the kSZ effect out to z ~ 1 using 
kSZ tomography This therefore provides a potentially 
powerful means of constraining our kSZ power spectrum 
models. 
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